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ABSTRACT 


It is shown that the plasma in Jupiter's ionosphere is col- 
lisionless above a certain level. In the outer magnetosphere where 
the rotational force dominates the gravitational force the collision- 
less plasma has a beam like distribution and gives rise to a two 
stream instability. ThlB leads to trapping of plasma in the centri- 
fugally dominated region of the magnetosphere. Plasma is lost via 
recombination. Equilibrium trapped particle densities are calculated 
by requiring a balance between trapping by wave-particle interaction 
and loss by recombination. The results are compared with recent 
observations from Pioneer 10. We indicate that the observations 
require an unexplained ion-heating mechanism. Some consequences of 
the model are discussed. 


1. INTRODUCTION 


A knowledge of the distribution of low energy plasma in the 
Jovian magnetosphere is of fundamental importance for our under- 
standing of the magnetic field structure and hence energetic particle 
measurements made by Pioneer 10 and 11. In particular, Smith et 
al. [1974] have shown that beyond 10 Rj the magnetic field differs 
markedly from that of a (tilted and displaced) dipole. They infer 
from their measurements the existence of a current sheet which 
requires the existence of considerable plasma densities. Frank t al. 
[1975] have reported measurements of low energy (E ~ 100 - 400 eV) 
protons which indicate relatively large densities of plasma between 
3 R_ and 15 R T . In section 2 we will review some of the published 
models for the distribution of low energy plasma in the Jovian mag- 
netosphere. We will point out that none of these are completely 
consistent with the results of Frank et al. Frank et al. 
observe protons with large pitch angles which must be trapped in 
the outer magnetosphere. Thus the model of e.g. Melrose [1967 ] which 
does not include trapped particles, is inconsistent with observations. 
Other authors ha.Va claimed that Coulomb scattering of photoelectrons, 
released f r -m Jovian ionosphere, will provide a sufficiently strong 
trapping mechanism. We believe, however, that this mechanism requires 
unreasonably high densities to be effective. And furthermore we 
will show in section 2 that the distribution function in the magnetosphere 


k 


i 

* 

of the particles released, from the ionosphere is beam like and unstable 
against a Landau resonant mode. Wave particle interaction due to this 
instability is a much more effective trapping mechanism than ordinary 
Coulomb collisions. 

Using quasi-linear theory we calculate the diffusion in 
velocity space arising from this instability. We will argue that the 
conditions of strong diffusion prevail in the Jovian magnetosphere. 

We then determine the total density in the Jovian magnetosphere 
by requiring that the loss of trapped particles (due to 
recombination) is balanced by the diffusion in velocity space. This 
is done in section 3. 

The resulting distribution function of trapped particles has 
a small perpendicular temperature but a rather large parallel tempera- 
ture, where pei-pendicular and parallel refer to directions relative 
to the magnetic field. This is clearly at variance with the results 
of Frank et al. who report rather large perpendicular temperatures 
(no measurement of parallel temperatures is available as yet). In 
section 4 we suggest two possible reasons for this. The trapped 
particles will, of course, undergo collisions which will tend to reduce 
the temperature anisotropy. The time scale for collisions is several 
hours and this process may be effective. Another suggestion follows 
the work by Kern [1962] who has shown that a sudden compression of 
the earth magnetosphere by the solar wind and the subsequent propagation 
of a compression pulse through the magnetosphere can lead to a large 
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increase of perpendicular energy. It is well known by now thst the 
Jovian magnetosphere is extremely variable in size and it may be 
thought that the intermittent compressions of the Jovian magneto- 
sphere result in an enhanced perpendicular energy density of the 
protons. This process is investigated in section 4 . Unfortunately 
both processes cannot account for the large energies observed inside 
the orbit of the moon lo. Finally, a brief discussion of some of the 
consequences of the proposed model plasma distribution will conclude 
the paper. To simplify the treatment we will neglect the tilt of the 
magnetic dipole with respect to the rotational axis and treat only the 
case of the dipole axis being aligned with the rotational axis. 
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2. DISTRIBUTION FUNCTIONS IN THE MAGNETOSPHERE 


In a frame which corotates with the planet particles experience 
a gravitational and centrifugal force. These forces may be derived 
from, a potential jd (geopotential) which can be written as 




2 

cos 



( 1 ) 


where GM T = 1.267 X 10 2 ^ cm^/s" 


[Anderson et al. 1974 ] 


R c = 0.958 Rj = 67,907 km 


[Michel and Sturrock 1974] 


n_ a 1.76 x 10“ 4 S" 1 

u 


[see e.g. Gledhill 1967, Piddington 1967, Melrose 1967, loarmidis and 
Brice, 1971, Michel and Sturrock 1974], 

If the plasma were in statistical equilibrium everywhere inside 
Jupiter’s magnetosphere the density would be distributed as 

n = n Q exp (-njftg/SKX) , 
where m = m. + m 

X G / ! 


( 2 ) 
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T = 


T + T, 
e i 


are the average mass and. temperature of the plasma. 

It was pointed out hy many authors that this equation predicts 
an increasing plasma density in the equatorial plane (\ = 0°) beyond 
r = 2.3 R, and it was concluded by everyone that such a distribution 
is unphysical. One must look for some physical constraints which 
limit the density in such a way that it decreases with distance in 
some regular manner. 

One group of authors [Gledhill, Piddington, Michel and Sturrock] 

limit the plasma density be requiring that the corotational energy 
1 2 2 

density NmG r ) not exceed the magnetic field energy density 

(B /8rr). For a magnetic dipole field this leads to a l/r decrease 

of the plasma density. The critical density in this model even at 

4 -3 

r = 10 R. is larger than 10 cm . It is questionable whether the 
J 

density ever reaches such a value at 10 Rj. If the density were as 
high as this one should have observed considerable distortions of the 
magnetic field inside 10 R_. This was, however, not the case [Smith 
et al. 1974]. At larger distances the field is not a dipole field 
any more [Smith et al. 1974] and without further knowledge of the 
magnetic field structure the above criterion cannot be used to predict 
the plasma density. 
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Other authors [Melrose, Ioannidis and Brice] drop the assumption 

of statistical equilibrium. And, indeed, Michel and Sturrock point out 

that the mean free path for collisions may be larger than the scale 

height even in the dense parts of the ionosphere. Thus we are forced 

to deal with a collisionless plasma and calculate the density variation 

from Liouville's equation. In particular the Jovian ionosphere will 

be essentially collisionless above a level hg where the density has 

decreased to a critical value N„ so that the mean free path X is 

b 

larger than the scale height H. For a hydrogen plasma this yields: 


H + H + 


11.4 T^/ 2 fkT L / 2 kT 

NgtoA H = mj4 G [h s ) 


(3) 


This reduces to 


< 1.5 T 

O 


oo 


if T is measured in [°K]. 

The height hg above the level of maximum density 
determined by the equation: 



is then 


« 




(5ss*l 
\ N s J * 


(?) 



9 


At this level we assum?* an isotropic Maxwellian distribution function. 
The density above this level can be calculated from 

n i - j 1 fVj ( V T n > dv i % ' (6) 

where f^ is the distribution function of the jth species. We require 
conservation of total energy K 




+ nyk + <1/ - K 


( 7 ) 


(where f> is an electrostatic potential) and the 1. adiabatic invariant 


il 



( 8 ) 


For the stationary case the density is then given as 


n. = 2tt 


]Jv: 
JJ X 


y dK du 



f . (K,(x ) — 1 dK dd 

v " G “ ^ 


( 9 ) 


In Figure la we show along a magnetic field line. There are 3 

S M 0 

important levels where the geopotentials are respectively, 

jl^ is the maximum value of the geopotential along the field lines . 
fig is the geopotential in the equatorial plane and is the geo- 
potential at the level h g . It is obvious that only particles with 
sufficient energy will be able to cross the level Thus to 
calculate the density n^ the range of integration for K is from 
n i.fi r + Z.fi to infinity, and the range for p, is from 0 to K - m.0_ - q 
And at the level Mwe have 

n JM = "js {jb($ ■ 41 + - ^JJ ' tl0) 

-K/kT 

where we have assumed fj(K;d) « e ** . Requiring charge neutrality 
at this level determines the electrostatic potential difference 
^ . The density can then be written as 


-*[$ ' ^J/2kT 

N = N c e 1 G GJ . {11 

M o 

The distribution function at this level is still isotropic and is 
shown in Figure lb. 

The density at the equator was calculated by Melrose [19^7] 


using essentially the same method. One obtains 
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jp 3 (1 - erf x^) exp x^ - (l - (l - erf Xg) exp x 


'M 



4 m (i - W* 1 x i • (ia) 


^ is the magnetic field, strength at the level M and B Q is the mag- 
netic field at the equator. The density Nq at Io's orbit (L = 6) is 
shown in Figure 2 for different temperatures T assuming an ionospheric 
plasma density of 5 x 10 cm” . Beyond lo's orbit Eq. (12) predicts 
a decrease of Nq/N^ * Clearly except for very large temperatures 

the magnetospheric densities predicted by Eq. (12) are too small to 
account for the observations. 

The distribution function at the equator consists of two beams 
moving from one hemisphere to the other and is shown in Figure lb. It 
is quite obvious that there will be a gap in the distribution function 
if 0^ - 0^ > 0 , which will be the case along field lines with 
L ^ 1.9739 [see Michel and Sturrock 197*0. The well known Penrose 
criterion predicts an electrostatic instability for this distribution 
function independent of the value of the temperature. (Any distribution 
function with a gap is unstable.) This instability can, however, be 
quenched by the presence of some plasma particles with velocities within 
this cap. These particles will not have a large enough energy or mag- 
netic moment to move back over the geopotential hill into the ionosphere. 
Thus any particle in the gap will be trapped in the outer magnetosphere. 



12 


Contrary to Melrose's prediction that there are no trapped particles 
we believe that there must be trapped particles on field lines with 
L ^ 1.9759. And, indeed, if there were no trapped particles the plasma 
detector onboard Pioneer 10 would have observed no plasma, because 
the detector is arranged in such a way that it can only measure 
particles with large pitch angles. Without trapped particles there 
would be none of these in the outer regions of the magnetosphere. 

Frank et al. [19753 have, however, reported considerable number 
densities of ions £-.1 magnetic field lines with L ranging from 
2 to 10. 

Michel and Sturrock [197^3 suggest that Coulomb scattering of 
photoelectrons is sufficient to populate the outer magnetosphere and 
thus provide for trapped particles. But Coulomb scattering of 
photoelectrons is a very ineffective process unless the density of 
the scatterers is large. Michel and Sturrock derive a criterion 
for trapping of photoelectrons by Coulomb collisions; namely that the 
mean free path for scattering is smaller than the length of a field 
line between the two positions where = 0^. For 10 eV photoelectrons 
this requires an average density (n,) along the field line between the 
two geopotential maxima 


<n c > 


1.5 X 10' 


cm 


It appears doubtful that photoelectric densities of this magnitude 
exist everywhere along field lines with L < 10. 
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Swartz et al. [1975] have recently determined the flux of 

photoelectrons escaping from the Jovian ionosphere. They find an 

average value for the total escape flux with energies from 3 eV to 
7 -2 -1 

100 eV as 10' cm s . The average energy of these particles is of 

-1 -3 

the order of 10 eV. Thus their average density is only about 10 cm 
which is much smaller than the critical density for trapping. In 
other words: the photo_lectrons will move from one hemisphere to 

the other without making a collision. The characteristic time for 

Q 

collisions is 10 s. But this does not mean that there is, on the 

8 8 

average, 1 collision per 10 s and hence 1 particle trapped in 10 s. 

If this were so one could rightly argue that trapping by collision is 
a cumulative process and the density of trapped particles becomes 
large after several collision times. The point is that the photo- 
electrons leave the trapping region (the region between the two 
geopotential maxima) before they can even make one collision. The 
flux lost by collisions along a field lines is 

F « F q (1 - e' X / X ) , (15) 

7 -P -1 

where F q ~ 10' cm s , x is a length along the field line and \ is 
the mean free path. The average total flux lost by collision is the 
flux of particles into the trapping region 
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r tr ' § £ F ** ” 8 ''of 1 + I - *>] * F c £ - M 

for small values of *£/X. For lo's flux tube (L = 6) 

X = 14.6 Rj = 1.2 X 10 11 cm C X . 

One may argue that although this is a very small flux it will 
eventually lead to a considerable density because the trapped par- 
ticles are lost only be recombination which is a very slow process. 
If recombination is the only loss mechanism, we can estimate the 

average density of trapped particles. The recombination loss per 

o -12 3 -1 

unit volume is a^ c If where w 5 X 10 cm s for a hydrogen 
plasm. The total l.oss due to recombination in the entire flux tube 
of area A is 


I A d^, N 2 dx = x a EC <N > 2 <A> 


(15) 


Thus we find 



F o 1 F o 1 <B) 

X a^, <70 X b m 


(16) 
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The mean free path for collisions which lead to trapping was given 
by Michel and Sturrock [1974] as 


X 


ON 



„ 2 
N a cos 
o 



(17) 


where a is the angle by which the particles have to be scattered in 
order to be trapped. This angle is the loss cone angle a which is 
given as 


sin a = b/b^ f 

ff = 2.49 X 10 cm 2 and p = v/c is the relative velocity of the two 
scatterers, i.e., of 2 particles coming from the southern and northern 
hemisphere respectively. Taking average quantities along a field line 
one finds the average density is 


F j?' 

(N) = 5 x 10 ^ -r cos 2 a = 2 x 10 ^ cos 2 a, [cm ^] (l8) 


E 


-2 -1 

where E is measured in eV and F in cm s 

o 

Using the values of Swartz et al. [1975] for the total flux 

escaping across the potential barrier, J* j(E) dE where E^ is the 

T1 



energy required to drag a proton across the harrier [see Xoannidis 
and Brice 19713/ we find that the equilibrium density is 0,5 cm"*^ 
at L = c 1 and 2 cm - ^ at L = 10. These are smaller than the values 
given by Xoannidis and Brice for two reasons. The flux of Swartz 
et al. is an order of magnitude smaller than the fluxes used by 
Ioannidis and Brice. Icannidls and Brice also assume that all photo- 
electrons will be trapped by collisions whereas we have shown that 

-6 

only a small fraction (typically of the order of 10“ ) will be trapped 
by collisions. 

It seems that Coulomb scattering cannot provide for an efficient 
trapping mechanism which is apparently required by the observations 
of Frank et al. It would only be an effective mechanism if X = T. 

We do not believe that this is the case inside 10 Rj. Instead, we 
believe tha. 'the essential physics Is contained in the instability 
created by the beams. This was already suggested by Ioannidis and 
Brice as well as by Goertz [1974] and Michel and Sturrock [1974]. 
However, no analysis of the effect of the instability on plasma 
densities has been published so far. 
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3. TRAPPING BY WAVE-PARTICLE INTERACTION 

The distribution function at the equator (as shown in Figure 
lb) is an inverted loss -cone distribution. It is known that this 
kind of distribution function is unstable [see e.g. Stringer, 1964]. 
Since the energy density of the beam particles is much less than the 
magnetic energy density we can treat the magnetic field as very 
large and only Landau resonances will be excited. The electric 
fields of the unstable waves are potential and there is no pertur- 
bation of the magnetic field. Then the wave-particle Interaction 
will cause only one-dimensional diffusion in velocity space along the 

characteristics v = const. [Kennel and Engelmarm 1966]. 

II 

Numerical plasma simulations have shown that for an infinite 
plasma the oscillation energy attains a value of about 5-10$ of the 
beam energy. In the case of Jupiter we do not, however, have an 
Infinite plasma but a finite one. Furthermore, we have constant 
injection of beam particles into the centrifugally dominated region 
of the Jovian magnetosphere. Tsytovitch [ 1970 ] has shown that under 
these circumstances the oscillation energy can reach very high levels 
through a pile-up effect- The reason for this is the following: 

Beam particles which enter the centrifugally dominated region 
tfS - > °> are not only acted upon by the waves they themselves 

generate but also by those excited earlier by other beam particles. 
Tills is so because the oscillation energy is convected along the field 
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lines with the group velocity dw/dk w which is smaller than the 
beam velocity. Newly arriving beam particles generate waves from an 
already high level created by the previous particles. This leads to 
e much higher level of oscillation energy. A calculation of this 
(level of oscillation energy is beyond the scope of this paper. 

Instead, we will argue that there exists a stationary level of electro- 
static oscillation in the centrifugally dominated regions of the 
magnetosphere. Beam particles and trapped particles will interact 
with this random force in 3uch a way that the loss rate due to recom- 
bination is balanced by the diffusion in velocity space. 

The analysis of the interaction proceeds in the same way as 
in the usual deviation of the quasilinear theory. We will deal with 
the one-dimensional case only. Then 

SF 

ds e" ikVS <^(t) E^t - s)> -J- ( 19 ) 

where 

f d ■ I p 3 \ d \ • (£0) 

f.(v ) is shown in Figure 5* 

2 II 
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For stationary oscillations 


(Ej^t) E^t - s)> » <^( 0 ) E^-s)) 


( 21 ) 


For simplicity we may assume an exponential correlation function 


(E^t + t) E^t)) = | E^. | 2 exp (-io^T - a k | t j ) 


( 22 ) 


with > 0 and a^ oj^. Then we obtain the diffusion equation 




|Ej ; 




m 2 ^ ) i(kv - u^) + a fc dv 


(23) 


We now replace the discrete spectrum by a continuous spectrum and 
consider the real part only 


!!i = ± 

dt dv 


%-/ hj* •<* - ^> «k g 
L d 


_d_ fD* Sf J 
dv l v dvl 


D* = I h 


m* \ 

J v 


(24) 
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It seems plausible that the spectrum of the oscillations is flat at 
least within a certain range of frequencies. Here we will use 
D* ** const. The total distribution function f^ consists of trapped 
particles fj and beam particles f^ supplied by the ionosphere. First 
we deal with the trapped particles. For a stationary state the loss 
afVa-t can be approximated by f^/T which assumes a recombination 
coefficient independent of velocity. Then a solution of Eq. (24) 
which is finite at the origin is 

r/(D* T^,) 1 / 3 , (25) 

where Is an Ayry function. According to the arguments above, D* 
is large and the argument of Bjf is small. 



f 5 ra0B i 


This indicates the formation of a plateau in the velocity distribution. 
This in in accordance with the general conc3usions of Kennel and 
Engelmann [19 66 ], The density of trapped particles is then 
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For large values of D* T this reduces to 


n. 



(28) 


Wj is the maximum velocity of the trapped particles 


no m *1* „ 2 kT . R, - B 

I («?) = $ - ta * (1 4 " 0) + - Z - 1 ~ 


2 '"J' ''G ''G m 


m. B 
J 


£n 


Pm" 


B 


? 2 kT . B m - 

\C. + - 3 . — p /n 

jo m^ B 


B 


B, 


'M 


B M “ 


( 29 ) 


For the beam particles we assume an ionospheric source of the form 


_b_ 

dv 



for lv I > W. and 
’ n 1 
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T =d ’ 2 “>3 • 

The life time t of the beam particles can be calculated as the ratio 
of the total beam density to the diffusion loss 



The mean free path of a beam particle is of the order of tv. or the 

J 

characteristic wavelength of the oscillations times the ratio of 

thermal energy to oscillation energy. As a rough estimate we take 

the wavelength as the beam velocity divided by the beam plasma 

frequency or 10 -1+ V^/F 1 / 2 . For F » IQ 6 cm” 2 s" 1 and V b » 10 8 cm 

this becomes 10^ cm and the mean free path is smaller than the length 

of a field line even if vhe energy of the oscillations is only l/lOO 

of the thermal energy. Under these conditions nearly all the beam 

is 

particles will be trapped. Then we can calculate the density by 
equating the loss due to recombination to the total flux provided by 
the ionosphere and not to only F (X/j) as in the case of Coulomb's 
collisions . 


M 


vm*- 4 .'; °kc< o \i> A(s)d8 


(32) 



where s = 0 In the equatorial plane. We anticipate that the trapped 

particle density is much larger than the beam density. Then charge 

o o 

neutrality requires that W“~ - W7 combining Eqs. (29) and (32) we 

6 1 

finally obtain for the trapped particle density 




( 33 ) 


where 



For the temperature of the beam particles we take the characteristic 

6 -2 -1 

energy of the escaping photoelectrons. For F^ = 10 cm s , 

— 1*2 

Q^ c = 5 X 10 and a surface field of 4 Gauss the trapped particle 
density in the equatorial plane is shown in Figure 4 . For L < 4 
the density must level off because the basic assumption that all the 
flux will be trapped is not valid anymore. This effect is difficult 
to calculate because it requires a knowledge of the diffusion constant 
D* and hence the spectrum of the electrostatic wave noise. Inside 
L = 2 no trapping occurs and the density will fall to values predicted 
by Eq. (12) which are very small. The solid curve in Figure 4 repre- 
sents the values obtained by Eq. (33 ), whereas the dotted curve is an 
indication of the effect of reducing the flux. The dashed line in 
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Figure 4 represents the density limit obtained by Gledhill [196? ] 
and Michel and Sturrock. [1974] for a surface field of 4 G. At h sn 30 
the plasma density will reach this critical value and one would 
expect serious distortions of the magnetic field beyond that distance . 
It is interesting to note that Goertz et al. [1974] find that the 
magnetic field observed by Pioneer 10 can be represented by the dipole 
field plus a perturbation field due to a current sheet which starts 
at 50 Rj. However, even at distance?, smaller than 30 Rj the drift 
motion of the plasma particle should provide for a current density 
and hence a distortion of the field. Then our method of calculating 
N*. which depends on the assumption of a dipole field, is not valid 
anymore. Thus the densities outside say 15 R T can not be considered 

U 

as very realistic. 


4 . COMPARISON WITH OBSERVATIONS 


Figure 5 shows equidenslty eontour plots inside 15 R_ calcu- 

U 

lated from « Cw^. The insert compares the values from the model 
with those inferred by Frank et al. from their observations. The 
agreement is quite good considering the assumptions made in the model 
and the uncertainties of the observations. Figure 5 displays clearly 
the disk- like shape of the plasma distribution first suggested by Gledhill 


[ 1967 ] which, however, he deriveu from somewhat different assumptions. 
The relatively dense plasma will cause a stretching of the field 
lines beyond 15 Rj which has the effect of making the disk thickness 
smaller because of the dependence of Cw, on B(s) and radial distance. 

U 

It should be obvious from the previous section that velocity 
space diffusion due to Landau resonances will leave the perpendicular 
energy of the particles unaffected. Thus we expect the trapped par- 
ticle distribution tc have a perpendicular temperature comparable to 
the ionospheric temperature. The parallel "temperature" should be 
much larger. Frank et al. have measured protons with characteristic 
energies of 100 eV. At present it is impossible to tell from their 
data whether the particle distributions are Isotropic or not. 

Assuming a Maxwellian distribution for the perpendicular velocities, 
we predict an anisotropic directional flux of the form 
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j(E,a) « AE e" E sin a / kT (1 + bE cos 2 a) (55) 

which, however, ia a strongly varying function of pitch angle (X only 

for small values of a. It is doubtful that an effective temperature 

of 100 eV could be explained by this dependence on a, particularly as 

the minimum value of a at which Frank et al. measured proton fluxes 

is about 4o°. Thus their characteristic energies are representative 

of the perpendicular energy and little can be said about the parallel 

energy. It should be noted that the spin modulation of the fluxes 

mentioned by Frank et al. may be related to the anisotropy described 

by Eq,. (56). It seems, however, difficult to separate this effect 

from the spin modulation due to penetrating high energy particles . 

We find it more likely that the relatively large energies reported 

by Frank et al. are indicative of a genuine heating of the proton's 

perpendicular energy. There are two possibilities which one might 

invoke to explain the large observed perpendicular energies . 

Trapped particles remain in the outer regions for a long time 

and are eventually lost only by recombination. The time scale for 

recombinations is ~ 2(a k)"\ For Ct , ~ 5 X 10 -12 cm^ s” 1 

He N e 

10 

this is of the order of 10 s at L = 6. The time scale for Coulomb 

collisions is only several hours. Clearly the trapped particles will 

undergo many collisions before they are lost by recombination. These 

collisions will eventually lead to an isotropic distribution with a 

1 2 

thermal energy of about 1/5 of the parallel energy w . Figure t 
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1 2 

shows - in the equatorial plane as a function of distance from 

the planet. We assume again that T is of the order of the energy of 
the escaping photoelectrons. These energies are generally of the 
right order of magnitude. But they are representative of the parallel 
energy. For an isotropic distribution the energies would be smaller 
than the observed ones by at least a factor of 6. The energy of the 
trapped electrons \ m w 3 is smaller than this by a factor of m /m 
if charge neutrality is valid. 

A second possibility of increasing the perpendicular energy 
may be a mechanism which was suggested by Kern [1962] for the earth's 
magnetosphere. He shows that repeated hydromagnetic shock waves 
accelerate particles trapped in the magnetosphere. For perpendicular 
shocks only the perpendicular energy is changed. He considers the 
shock waves as being generated by a sudden compression of the front 
side of the magnetosphere by increases in the solar wind pressure. 

At some distance in the magnetosphere the solar wind compresses the 
magnetosphere faster than the local hydromagnetic wave velocity and 
the compression pulse propagates as a shock wave. If the rise time 
of the shock is smaller than an ion gyroperiod ions will not conserve 
their adiabatic moment, but the magnetic moment of a trapped particle 
is increased as it is swept up by the shock. The subsequent slow 
relaxation of the compression conserves the adiabatic moment and after 
the magnetic field returns to its initial value the perpendicular 
energy (E^ = UB) wil\ oe larger than before th^ shock event. The energy 
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gain of an ion by a single shock is not large enough to explain 
Frank's result. But if particles are subject to a number of shocks 
it is conceivable that their perpendicular energies increase con- 
siderably. Kern finds that proton energies of up to several keV can be 
obtained by this process in the earth's magnetosphere. 

If the time between successive shocks is large compared with 
the time for collisions we would not expect this process to be 
important. We do, however, know by now that the Jovian magnetosphere 
is extremely variable in size and that the solar wind does indeed 
push the magnetosphere in. Several events have been witnessed by 
both Pioneer 10 and 11 [see e.g. Smith et al. 1975, Wolfe et al, 1974], 
Kern finds that the transverse kinetic energy .tensity e n in the 
undisturbed medium after n shock events is increased above the trans- 
verse energy density in the undisturbed medium by 


e 

n 


e 

o 



N 1/5 

Kl 


2n 


(56) 


where U Q is the velocity of the compression pulse and is a local 
Alfven velocity. For a shock we must have U Q > C^. Since this 
process leaves the density unchanged we would expect the perpendicular 
energy to increase as 


fro* 1 

r°ii 

1/5" 

ki + 1 

l u oJ 



£n 


E = E 2 
i xo 
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beyond the distance where U Q -> C^. Since the Alfven velocity 
(C^ « B/ _y4nm i N ) decreases with distance this process may conceivably 
lead to an increasing perpendicular energy. However, a typical 
Alfven velocity inside r = 10 Rj is 2 x 10 cm/s, and one would 
require very fast compressions for this mechanism to work. There 
is little hope that this mechanism can account for the 100 eV protons 
observed inside lo's orbit. It may be active in the distant parts of 
the magnetosphere and in part responsible for the energetic protons 
required for the observed perturbation of the magnetic field. 

Since there is very little doubt that Frank's energies refer 
to perpendicular energies the above arguments show that there must be 
a mechanism active in the Jovian magnetosphere which heats the ions. 
If, as Frank et al. assume, the observed ions come from the Jovian 
ionosphere and have indeed a temper attire of 100 eV everywhere along 
the field lines, there must be a very large density of extremely hot 
ions in the ionosphere. Then the geopotential barrier is relatively 
unimportant, the plasma is collision dominated everywhere and the 
treatment above is irrelevant. But then one encounters the usual 
difficulties of a model in which the plasma is in statistical equi- 
librium (see section 2). On the other hand, the assumption of an 
outside source (say the solar wind) and radial diffusion would not be 
able to explain the fact that the energy decreases with decreasing 
distance unless severe energy losses were involved. Furthermore, the 
solar wind could not supply the large densities observed. We believe 
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that the energization of the ions must take place along field lines. 
The following is a very speculative attempt to describe such a 
mechanism. 

In the previous section we have assumed charge neutrality 
everywhere along a field line. This may not be a valid assumption. 

It has been speculated that particle acceleration along field lines 
in the earth auroral ionosphere is due to potential double layers 
[see e.g. Block 1972 and references therin]. Potential double 
layers are small regions where charge neutrality is not valid and 
strong electric fields exist. 

These double layers, or electrostatic shocks, presumably form 
when a mechanism exists through which an appropriate distribution 
function for trapped particles (i.e., those which are reflected by the 
double layer) can be maintained. Recent computer simulations [Goertz 
and Joyce 1975] have shown that stable double layers can be generated 
under certain circumstances. Knorr and Goertz [197^-] have described 
double layers in the framework of BGK solutions of the Vlasov equation. 
They find that strong double layers require distribution functions 
for trapped particles which have a minimum at v = 0, precisely those 
cue would obtain in the Jovian magnetosphere. The potential drop 
across a double layer may well obtain many KT/e as in the numerical 
simulations. Two double layers, one in the northern hemisphere and 
one in the southern hemisphere, would accelerate the ions and at 
the same time prevent them from escaping along field lines into the 


opposite hemisphere. Figure 7 is a sketch of the possible location 
of these structures. Although this model may explain the large 
observed energies much more work is needed in order to assess the 
possibilities of the various mechanisms mentioned. However, this 
seems impossible without any further information about the pitch 
angle distribution of the trapped protons. 
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5. DISCUSSION 

In the previous sections we have shown that the Jovian 
ionosphere can provide enough plasma to explain the observations of 
Frank et al. We also showed that an ion-heating mechanism seems to 
be required. Note that as long as a strong diffusion is valid the 
trapped particle density is independent of the characteristic 
temperature of the trapped particles . In this section we want to 
discuss two consequences of our model, one relating to the densities 
only and the other relating to the large characteristic energies of 
the trapped particles. 

It is known [Kennel and Petschek 19 66 ] that the stabley 

trapped high energy particle fluxes arc. determined by the intensity 

of whistler mode noise. Particles above a certain critical energy, 

2 

which scales like B /8 ttN, are pitch angle scattered and lost more 
rapidly to the ionosphere than particles with smaller energies . 

Figure 8 shows how B /8 ttN varies with distance (assuming a dipole 
field of 4 G surface strength). These values are not unlike the 
values at which Baker and Van Allen [1975 3 and Filius et al. [1974-] 
observe a change in the energy spectrum of the electrons. A detailed 
analysis of the energy spectra is being undertaken at present and 
the results will be published elsewhere. 
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Because the characteristic energy of the trapped particles 
is rather large , certain ideas about the distribution of plasma in 
the magnetosphere may not be valid. In particular, Gledhill [1967] 
Goertz [1972] and Hill et al. [197*0 have predicted that the thermal 
plasma should be confined to a centrifugal symmetry surface which is 
inclined with respect to the magnetic equatorial plane but which is 
not parallel to the rotational equatorial plane either. This is due 
to the fact that the centrifugal force tends to push the plasma to 
the most distant point (from the rotational axis) along a field line. 
This argument is, however, only valid if the centrifugal force is 
much larger than the magnetic force, i.e., 


1 n 2 . 2 2 . . 1 
2 n ?(r cos » -g 


1 2 2 
2 V i + V || 
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This inequality is true for a low-temperature plasma. But in the 
Jovian magnetosphere the energy of the plasma is at least comparable, 
if not larger, than th^ corotational energy mil r /2. Thus, the 
inequality is not fulfilled in the Jovian magnetosphere. Indeed, it 
seems much more likely that the inequality i3 actually reversed. 

Then the particles are confined to the magnetic equatorial plane. 
Without a knowledge of the pitch angle distribution of the trapped 
particles we cannot decide whether the inequality is reversed or not. 
There is some experimental evidence, derived from Pioneer 10 magnetic 
field data [Goertz et al. 197*0, which indicates that the thermal 
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plasma is confined to the magnetic equatorial plane. Clearly, more 
information on the plasma distribution and temperatures is needed 
in order to decide these questions and assess the validity of the 
model presented. It is our sincere hope that ft proper plasma experi- 
ment will be Included on future Jupiter probes ii 




35 


ACKNOWLEDGMENTS 


It is a pleasure to acknowledge helpful discussions with 
Professors J. A. Van Allen and L. A. Frank, as well as with Dr. 

K. Ackerson. Some of the ideas expressed in this paper originate 
from discussions with Professor J. A. Gledhill. 

This work was supported in part by NSF Grant GA -31676 and 
NASA contracts NAS2-5603 and NAS2-^,553. 


36 


REFERENCES 

Anderson, J. D,, G, W. Null and S. K. Wong, Gravity Results from 
Pioneer 10 Doppler Data, 79; 3661 (197*0 * 

Baker, D. N. and J. A. Van Allen, Energetic Electrons in the Jovian 
Magnetosphere, Univ. of Iowa Research Report 75-13 (1975). 

Block, L. P., Potential Double Layers in the Ionosphere, Cosmic 
Elect rodyn . , 3, 34 9 (1972). 

Filius, R. W. and C. E. Mcllwain, Measurementr of the Jovian Radiation 
Belts, J. Geophys. Res., 7£, 3589 (1974). 

Frank, L. A., K. L. Ackerson, J. H. Wolfe and J. D. Mlhalov, 

Observations of Plasmas in the Jovian Magnetosphere, Univ. 
of Iowa Research Report 75-5 (1975). 

Gledhill, J. A., Jupiter's Magnetosphere, Nature , Lond., 214, 156 (1967). 

Goertz , C. K., Ph.D. Thesis, Rhodes University, Grahamstown, South 
Africa (1972). 

Goertz, C. K. , Jupiter's Ionosphere and Magnetosphere, Planet. Space 
Sci., 13, 95 (1974). 

Goertz. C. K. and G. Joyce, Numerical Simulation of the Plasma 
Double Layer, Astrophys. and Space Sci ., 32, 165 (1975). 

Goertz, C. K,, R. A. Randall and M. F, Thomsen, A Model for the 
Jovian Magnetosphere, EOS , 56, 1172 (SM50) ( 1974 ). 


Hill, T. W,, A. J. Dessler and F. C. Michel, Configuration of the 

Jovian Magnetosphere, Geophys . Res . fetters , 3 (1974). 

Ioannidls, G. and N. Brice, Plasma Densities in the Jovian Magnetosphere: 

Plasma Slingshot or Maxwell Demon?, Icarus , 14, 3fQ {1971). 

*■ 

Kennel, C. F. and F. Engelmann, Velocity Space Diffusion from W 
Plasma Turbulence in a Magnetic Field, Phys. of Fluids , 9, 

2377 (1966). 

Kennel, C. F. and H. E. Petschek, Limit on Stably Trapped Particle 
Fluxes, J, of Geophys . Res . , 71 , 1 (1966). 

Kern, J. W. , A Note on the Generation of the Main-Phase Ring Current 
of a Geomagnetic Storm, J. of Geophys. Res ., 67, 3?37 (1962). 
Knorr, G. and C. K. Goertz, Existence and Stability of Strong Potential 
Double Layers, Astrophys, and Space Sci ., 31, 209 (197*0 • 

Melrose, D. B., Rotational Effects on the Distribution of Thermal 
Plasma in the Magnetosphere of Jupiter, Planet. Space Sci . , 

15, 381 (1967). 

Michel, F. C. and P. A. Sturrock, Centrifugal Instability of the 

Jovian Magnetosphere and its Interaction with the Solar Wind, 
Planet. Space Sci ., 22, 1501 (1974). 

Piddington, J. H. / Jupiter 's Magnetosphere, Univ. of Iowa Preprint 
67-63. 

Smith, E. J., L. Davis. Jr., D. E. Jones, P. J. Coleman, Jr., D. S. 
Colburn, P. Dyal, C. P. Sonett and A. M. A. Frandsen, The 
Planetary Magnetic Field and Magnetosphere of Jupiter: Pioneer 

10, J. Geophys . Res . , 79, 3501 (197U). 


Smith, E. J., L, Davis, Jr., D. E. Jones, P. J. Coleman, Jr., P. 

Dyal, C. P, Sonett and A. M. A. Frandsen, Jupiter’s Magnetic 
Field, Magnetosphere and Interaction with the Solar Wind, 
to be published in Science (1975). 

Stringer, T. E., Electrostatic Instabilities in Current-Carrying 

and Counter-Streaming Plasmas, J. Nucl. Energy . l6, 067 (19^4) . 

Swartz, E. W. R. W. Reed and T, R. McDonough: Photoelectron Escape 

from the Ionosphere of Jupiter, J. of Geophys . Pe s., t 80 , 495 (1975). 

Wolfe, J. H., J. D. Mihalov, H. R. Collard, D, D. McKibbin, L. A. Frank 
and D. S. Intriligator, Pioneer 10 Observations of the Solar 
Wind Interaction with Jupiter, J. Geophys . Res . , 79, 3489 (197^)* 


39 


Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

i ; 

Figure 


FIGURE CAPTIONS 

la The variation of the geopotential along a field 
line. The curve is not drawn to scale. 

lb The distribution functions of particles from the 
ionosphere at the levels S > M and 0. 

ii 

2 Collisionleas densities in the equatorial plane at 
6 Rj for different ionospheric temperatures . 

3 The parallel distribution function f . « u r F v dv . 

d j x i 

4 Equatorial trapped particle densities. For a description 
of this figure refer to the text. 

3 Equidenslty contour plots. The contours are in intervals 
of 10 cm"*^. This figure does not take into account the 
tilt of the dipole with respect to the rotational axis. 

6 Characteristic energies of trapped particles. The dots 
indicate some "temperatures" observed by Frank et al, [19751* 

7 A sketch of the possible location and polarity of potential 
double layers in the Jovian magnetosphere. 

8 The magnetic energy per particle (B^/8 ttN) in the equatorial 
plane as a function of distance from Jupiter. 
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